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Abstract 

It is shown that unlike the perfect fluid case, anisotropic fluids 
(principal stresses unequal) may be geodesic, without this implying 
the vanishing of (spatial) pressure gradients. Then the condition of 
vanishing four acceleration is integrated in non-comoving coordinates. 
The resulting models are necessarily dynamic, and the mass function 
is expressed in terms of the fluid velocity as measured by a locally 
Minkowskian observer. An explicit example is worked out. 
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1 Introduction 



As is well known, the vanishing of four acceleration (geodesic condition) 
implies for a perfect fluid that pressure gradients vanish. In the case of 
spherical bounded ( non-dissipative) configuration, the vanishing of pressure 
at the boundary surface, implies in turn the vanishing of pressure everywhere 
within the distribution (dust). 

Indeed, for a perfect fluid the equation of motion reads 

(p + p)a a = h a »p tU (1) 

with 

h% = $Z- u a u, (2) 

a» = u v u% (3) 

where colon and semicolon denote partial and covariant derivatives, and as 
usual a M , p and p stand for the four acceleration, the four velocity, the 
energy density and the pressure respectively. 

^From the above it becomes evident that the geodesic condition implies 
the vanishing of pressure gradients. From purely physics considerations this 
conclusion is also obvious: the vanishing of four-acceleration means that only 
gravitational forces are acting on any fluid element, thereby implying that 
pressure gradients (the only hydrodynamic force in a perfect fluid) vanish. 
However in the case of anisotropic fluids, an additional force term appears 
besides the pressure gradient (see next section). Therefore it is in principle 
possible to have a fluid distribution, such that both terms cancel each other, 
leading to a geodesic fluid with non-vanishing pressure gradients. 

Since the original Lemaitre paper and particularly since the work of 
Bowers and Liang anisotropic fluids have attracted the attention of many 
researchers in relativity and relativistic astrophysics (see || and references 
therein), due to the conspicuous role played by local anisotropy of pressure 
in the structure and evolution of self-gravitating objects. It is the purpose 
of this work to present further models of anisotropic spheres, based on the 
geodesic condition. Besides the natural interest of such models in general rel- 
ativity, the presented models are interesting because they represent the gen- 
eralization of Tolman-Bondi [|J models to anisotropic fluids, in non-comoving 
coordinates. Incidentally, it is worth noticing that in the classical paper by 
Oppenheimer and Snyder on dust collapse [|], they start their study, using 
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the same kind of coordinates we use here, and then switch to comoving ones, 
in order to integrate the field equations. 

The plan of the paper is as follows. In Section 2 we define the conventions 
and give the field equations and expressions for the kinematic variables we 
shall need, in noncomoving coordinates. The geodesic condition is explicitly 
integrated in Section 3. In Section 4 we work out an example. Finally a 
discussion of results is presented in Section 5. 

2 Relevant Equations and Conventions 

We consider spherically symmetric distributions of collapsing anisotropic 
fluid, which we assume to evolve adiabatically (without dissipation), bounded 
by a spherical surface E. 

The line element is given in Schwarzschild-like coordinates by 



where v{t,r) and \(t, r) are functions of their arguments. We number the 

coordinates: x° = t; x 1 = r; x 2 = 6; x 3 = 0. 

The metric @ has to satisfy Einstein field equations 



ds 2 = e u dt 2 - e x dr 2 - r 2 (d6 2 + sin 2 Qd<f) , 



(4) 
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which in our case read 0: 
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where dots and primes stand for partial differentiation with respect to t and 
r, respectively. 

In order to give physical significance to the components we apply the 
Bondi approach ||. 

Thus, following Bondi, let us introduce purely locally Minkowski coordinates 
(r, x, y, z) (alternatively one may introduce a tetrad field associated to locally 
Minkowskian observers) . 

dr = e u ^ 2 dt ; dx = e x ^ 2 dr ; dy = rdO ; dz = rsinOdcf). 

Then, denoting the Minkowski components of the energy tensor by a bar, we 
have 



rpl) 

1 n 



J ) 



fl = Tl ; f | 



T, 



2 > 



T? = T, 



3 ) 



f 01 = e ~ {u+x) ' 2 T { 



01- 



Next, we suppose that when viewed by an observer moving relative to these 
coordinates with proper velocity u)(t, r) in the radial direction, the physical 
content of space consists of an anisotropic fluid of energy density p, radial 
pressure P r and tangential pressure P±. Thus, when viewed by this moving 
observer the covariant tensor in Minkowski coordinates is 



( P 




V o 





Pr 










P± o 

P, 



Then a Lorentz transformation readily shows that 

p + P r oo 2 



rpO 
1 



rpO 
J 



Tl = Tl 



l-LU 2 ' 

P r + puj 2 

l-UJ 2 



(10) 
(11) 
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T 2 = T, 
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(p + P r )uje 
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Note that the coordinate velocity in the (t, r, 9, 0) system, dr/dt, is related 
to ijj by 

^,r) = |e^)/ 2 . (14) 



Feeding back (flOHl3l) into we get the field equations in the form 

p + P r u 2 1/1 _ A 



1 — lu 2 8tt [ r 2 

Pr +Pt0 2 _ 1/1 

1 — uj 2 8n\r 2 
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A ^z/" + z/ 2 - AV + 2^—^) }, (17) 



(p + P r )uje^ +X ^ 2 A 



STrr" (18) 
At the outside of the fluid distribution, the spacetime is that of Schwarzschild, 
given by 



ds ? = ^ _ 2^ dt 2 _ A ^ 1 rfr 2 _ r 2 ^2 + sm 2^ 



(19) 



As is well known, in order to match smoothly the two metrics above on the 
boundary surface r = r%(t), we must require the continuity of the first and 
the second fundamental form across that surface. In our notation this implies 

2M 

= 1 , (20) 



2M 



(21) 



and 

[P r ] s = 0, (22) 
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Where, from now on, subscript £ indicates that the quantity is evaluated at 
the boundary surface E. 

Eqs. (^0]), ( PH ) and (|22|) are the necessary and sufficient conditions for a 
smooth matching of the two metrics (f|) and (|19|) on E. 
Next, let us write the energy momentum tensor in the form 



TflU = (P + P±) U^Uy - P±g^ U + (P r - P±) S^Sy 



with 



-v/2 



lu e 



-A/2 



(1 -uj 2 ) 1 ' 2 ' (l-UJ 2 ) 
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-v/2 



-A/2 



(1_ ^2)1/2' (1 _ ^2)1/2'°' y 



(23) 

(24) 
(25) 



where w' 1 denotes the four velocity of the fluid and s M is a radially di- 
rected space-like vector orthogonal to u^. Then the radial component of the 
conservation law 



T M = 0. 



(26) 



may be written as 



(-an?)' 



16n 



A 2 Az> s 



T\-Tl WW 2? -2* +— A + --- , (27) 



which in the static case becomes 



2 r 



(28) 



representing the generalization of the Tolman-Oppenheimer-Volkof equa- 
tion for anisotropic fluids ||. Thus, as mentioned before, local anisotropy 
introduces an extra term in this "force" equation, besides the usual pressure 
gradient term. 

Finally, for the two non-vanishing components of the four acceleration, 
we easily find 



a 



UJUJ 
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(30) 



3 Integrating the Geodesic Condition 



Let us now integrate the geodesic condition. First, observe that from the 
field equations (|II!|),(|I6"D and (|TB|), one obtains after simple manipulations 



{u-X)/2 (V , / 



(\' + i/) + (l + u 2 )X = 



(31) 



Next, it follows at once from (|29|) and (0) that, 

uja x = -a e (x -" )/2 (32) 
Therefore the vanishing four-acceleration condition amounts to 



UJUJ V 



coX Co 



- + e -'W — + 



0. 



(33) 



Then, replacing v' by its expression from (|3l"l) , into (|33|) , this last equation 
becomes 



Lue 



or, using ( fT4[) 
whose solution is 



1 — or 1 — ur 

4>dt + = 0, 



<p = ln(l — uj 2 ) + A = Constant. 
Finally, from the fact that co(t, 0) = we obtain, 



e- x = l-uj 2 . 



Introducing the mass function as usually, 



2m 



(34) 

(35) 
(36) 

(37) 
(38) 
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we have 

m = — (39) 

In all the above we have assumed uj 7^ 0, since from simple physical 
considerations we should not expect static solutions to exist. 

Indeed, if we assume staticity (u = 0) then the geodesic condition implies 
v' = 0, which in turn, using flTED and (^), leads to 

2m , . 

8vrP r = - — . (40) 

Then junction condition (p2j) would lead to = M = 0. 

There is however one possible case of static geodesic solution, which ap- 
pears if we relax the condition of continuity of the second fundamental form 
(implying the continuity of radial pressure) across the boundary surface, and 
assume the existence of a surface layer 0. 

In this specific case, it follows from (40), the geodesic condition and field 
equations (15) and (17) that 

p + P r + 2P ± = (41) 

implying that the active gravitational mass (Tolman,||) defined for any 

r < rs as, 

m T = 4tt f r 2 e^ +X)/2 (T ° - 2* - 2T 2 )dr (42) 
Jo 

vanishes inside the sphere. 

We shall not consider here these kind of solutions and accordingly all our 
models will be dynamic (u 7^ 0) and satisfy all junction conditions . 

Now, from (|HD evaluated at the boundary surface, and fl20"D, (pl|), we 
obtain 

uj-£ = ^— (43) 



On the other hand, from (39) evaluated at the boundary surface, we have 



2M , . 

^ = ±y— (44) 
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Where the + (— ) refers to the expansion (contraction) of the surface (from 
now on we shall only consider the contracting case). Feeding back ( f44|) into 
,, we get 



Then equating (fi5D and fl44|) we have 

2M\ 3/2 



2M\ 1/2 



This equation may be integrated to give 



t 



2M 



2tanh 



I2M _ 2 [1 + 6M/r s ] 
r s 3 (2M/r s ) 3/2 



(45) 



(46) 



(47) 



giving the evolution of the boundary surface. Unfortunately this last equation 
cannot be inverted (at least we were unable to do that) to obtain the explicit 
form r = r^(t). Accordingly we have also integrated (TO) numerically , in 
order to exhibit the evolution of ru, see figure(l). 

So far we have found all consequences derived from the geodesic condition 
which, obviously, are valid in the pure dust case as well as in the anisotropic 
case. In the next section we shall work out an explicit example by impossing 
an "equation of state" for the physical variables. 



4 A Model 

The purpose of this section is not to model any specific physical system, but 
just to illustrate the consequences derived from the geodesic condition. Thus, 
somehow inspired by the incompressible fluid model, let us assume 

T ° = f{t) (48) 

then from ( |3"9"| ) and the fact that 

m > = 47rr 2 T ° (49) 

one obtains 
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where 



/(*) 



3M 
47rr| 



(51) 



Observe from (|50|) that the evolution in this model is homologous. 
Next, introducing the dimensionless variables 



x = 



we have 



y = 



i- x - 



2M 



u = —- 



X 



Finally, from the field equations the following relations follow 
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(57) 



Then from (55) and (p6|) we obtain p + P r as function of x and y. 
Feeding back this expression into (0),this equation may be solved for u, 
which in turn allows to express all physical variables (p, P r and P±) in terms 
of x and y which are given by (|47|) or alternatively by the numerical solution 

of " 



5 Conclusions 

We have seen that the geodesic condition, which for anisotropic fluids is 
compatible with the presence of pressure gradients, can be integrated, giv- 
ing the explicit form of the evolution of the boundary surface. Resulting 
models may be regarded as generalizations of Tolman-Bondi solutions, to 
anisotropic fluids. In order to obtain the evolution of all physical variables 
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for different pieces of matter, additional information has to be given. In the 
model above we have assumed condition (f48|), which in turn leads to the 
homology condition (pUp. Parenthetically, this last condition is widely used 
by astrophysicists in their modelling of stellar structure and evolution ||. 

Figures (1) and (2) display the behaviour of the radius and the evolution 
of ojy, in the contracting case. As expected, as the boundary surface ap- 
proaches the horizon, its coordinate velocity (r^) stalls, whereas the velocity 
ws measured by the locally Minkowskian observer, tends to light velocity. 
Figure (3) shows the sensitivity of the pattern evolution with respect to the 
compactness of the initial configuration. As expected more compact config- 
urations collapse faster. The good behaviour of p and p r is easily deduced 
from (0)-fT7p. 
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Figure captions 

Figure 1. y — r s /2M as function of t/M for the inital value y(0) = 30. 

Figure 2. uo^, as function of t/M for the same initial data as in figure 
I. 

Figure 3. c^s as function of t/M for y(0) = 30, 29, 28, 27, 26 curves from 
rigth to left respectively. 
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Figure 1 . 




Figure 3 . 



